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Abstract 

Studies on 'strange' particle production have always occupied a very important space in the 
domain of Particle Physics. This was and is so, just because of some conjectures about specially 
abundant or excess production of 'strange' particles, at certain stages and under certain condi- 
tions arising out of what goes by the name of 'Standard' model in Particle Physics. With the 
help of Hagedornian power laws we have attempted to understand and interpret here the nature 
of the pT-spectra for the strange particle production in a few high energy nuclear collisions, some 
interesting ratio-behaviours and the characteristics of the nuclear modification factors that are 
measured in laboratory experiments. After obtaining and analysing the final results we do not 
confront any peculiarities or oddities or extraneous excesses in the properties of the relevant 
observables with no left-over problems or puzzles. The model(s) used by us work(s) quite well 
for explaining the measured data. 
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1 Introduction 

Studying the nature of particle production in proton-proton collisions is important and interesting 
in itself, as it might shed light on the basic mechanism for production of particles. Besides, it could 
also serve as a necessary benchmark for the physics developments in ultra-relativistic heavy ion 
collisions [TJ [2]. This is specially important at the large hadron collider (LHC) where the heavy ion 
programme had started by November 2010 delivering some preliminary results on some aspects of 
strange hadrons produced as the final product in high energy nuclear collisions and these strange 
secondaries are supposed to provide valuable insights into the properties of the 'controversial' system 
newly formed. One of the main motivation for measuring strange particles in heavy nucleus-nucleus 
reactions at LHC is the expectation that their production-rates for participating nucleon should 
be enhanced with respect to basic nucleon-nucleon interactions. Strangeness enhancement has 
consistently been proposed as one of the strong diagnostics for a Quark Gluon Plasma (QGP) 
state O 0]. The enhancement factor (E) is defined as rapidity-density of multiplicity (yield) par 
mean number of nucleon participants [< N par t >] in heavy ion collisions, divided by the respective 
value in p+p collisions. The requisite information about < N par t > etc are to be obtained from 
Abelevet al 0E]. 

As the strange hadrons are not at all present in the initial system(A), the question rises very 
sharply: how do they make their appearance in the final products. So there must be some specific 
reflections on the constituent pictures of the particles and specifically the nucleons. Besides, en- 
hancement of strangeness productions was/is one of the powerful diagnostics for the formation of 
quark gluon plasma (QGP) in relativistic heavy ion collisions and the colliders(RHIC).The observa- 
tions of an increase of strange baryon production relative to p+p collision in SPS data, confirmed 
later at the RHIC studies, has brought excitement in this area. Besides, the increase of p/ir ratio 
(B) in such collisions in the non-strange domain had its parallel in the strangeness sector with the 
observation of slow rise of the A/A; values. 

The organization of this work is as follows. In the section-2 we give an outline of the model to 
be applied. In the next section(section-3)we deliver the results by figures and tables with a short 



discussion on the results obtained. In the last chapter, we precisely point out the conclusions to be 
arrived at. 

2 The Background In Some Detail And The Working Formulae 

With gradual attainments of larger transverse momenta(j>T) of the secondaries in high-pr(hard) 
interactions, the problems of deviations from exponential nature of fits on invariant spectra began 
to crop up steadily. Gazdzicki and Gorenstein [7J observed rightly that for px > 2 GeV/c, the 
data sharply deviate from the exponential nature, for which Darriulatl8| proposed a power law 
distribution of certain forms for both p-p-spectra an d particle multiplicity. Indeed, for both pr- 
spectra and multiplicity such power law forms have become now the most dominant tools in dealing 
with the transverse momentum spectra of all hadrons. Gazdzicki and Gorenstein showed that the 
normalised multiplicities and (jtit ~ Pt) spectra of neutral mesons obey the mr-scaling which 
has had an approximately power law structure of the form ~ (m-p)™, where my is called transverse 
mass and is defined by vht=\/p\ + tt? . This scaling behaviour was analogous to that expected 
in statistical mechanics, the parameter n plays the role of temperature and any normalization 
constant to be used resembles the system volume. Thus the basic modification of the statistical 
approach needed to reproduce the experimental results on some hadron production process in 
p(p) + p interaction in the large m? = Pt domain is to change the shape of the distribution 
functions exp(— — ) had to be altered to the power law form as given by (— )~ n with some changed 
parameters, viz, a scale parameter A and an exponent n, both are assumed to be common for all 
hadrons. 

Let us now dwell, in brief, on the clues to the possible origin of power laws. One of the ba- 
sic features of the hadronuclear collisions is: irrespective of the initial state, agitations caused by 
the impinging projectile(be it a parton or particle/nucleon) generate system effects of producing 
avalanches of new kind of partons (called gluons) which form an open dissipative system. And these 
production processes are not at all gradual; rather they are very sudden, drastic and complex. 
And such complex properties and processes in nature do generally subscribe to the power-law be- 
haviours. In the recent times, it is being propounded consistently that the power law behaviours 



put into use here are "manifestations of the dynamics of complex systems whose striking feature is 
of showing universal laws characterized by exponents in scale invariant distributions that happen 
to be basically independent of the details in the microscopic dynamics" |9j. The avalanches caused 
by production of excessive number of some new variety of parton called 'gluons'(the process called 
'gluonisation') give rise to the jettiness of particle production and of cascadisation of the particle 
production processes leading to the fractality as is shown by Sarcevic|10j. These cascades are self- 
organizing , self-similar and do just have the fractal behaviour. Driven by the physical impacts of 
these well-established factors, in the high energy collision processes do crop up the several power- 
laws. And how such power laws do evolve from exponential origins or roots is now-a-days being 
taken care of by the induction of Tsallis entropy [UJ and a generalisation of Gibbs-Boltzmann 
statistics for long-range and multifractal processes. 

In what follows we are going to choose a specific form of power law which was previously applied 
by us in the case of hadron-nucleus collisions. With a view to accommodating some observed facts 
for strange particle production, it is tempting to try to fit the whole distribution for the inclusive 
PT-spectra with one single expression in the form of power law as was done by G. Arnison et al|12j 
and HagedornfTB"]. 
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where the letters and expressions have their contextual significance. This parametrization describe 
the data well over the entire range of pt- 

Indeed for px — > 0, oo we have, 
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Thus along with impressive fit, which now includes the large pt domain, the estimate of < pt > 
assumes with the help of expression(l) : 

>= [ q/{PT + q) n p T dpT = 2q 
/ q/(PT + q) n PTdpT n - 3 



So, in clearer terms, let us put the final working formulae here as follows with substitution of px 
(transverse momentum) as x in the power-law model [14} \TE\ respectively 

f(x)=A(l + x/q)- n (4) 

There is yet another very important observable called nuclear modification factor (NMF), denoted 
here by Rcp, which for production of any hadron is defined by [16] 

ff , s [(d 2 N/(27rdpTdp T dy)) c /N bin } Cenrral 
CP[PT) [(d*N/(2ird PT dpTdy)) c /N bi n] Peripheral ( ' 

3 Results 

In obtaining the results presented here, no serious statistical calculational procedure was adopted. 
The graphs are drawn more as fitological-cum-phenomenological exercises with mainly statistical 
errors in considerations. The experimental data do not provide, in the most cases, any systematic 
errors. Data points for the heavy, high strangeness-valued particles are too scarce; for which the 
number of degrees of freedom is too limited for many cases. The quality of fits to the data indicated 
in the tables by x 2 /ndf terms in the columns is attempted to be kept at a modestly satisfactory 
value (tending as nearly as possible to unity). And the figures are drawn by Wgnuplot, wherein 
there are some inbuilt statistical procedures and techniques. 

Quite observably, the results are depicted here in graphical plots. And the used values of the 
corresponding parameters for fits are shown in separate tables. The Fig. la and Fig. lb show the 
production of secondaries k°,k + ,k~,A and A in proton-proton collisions at y / s=200 GeV at the 
rapidity y < 0.5(Table-l).The figures in Fig. 2a and Fig. 2b depict the results for the H~,H + ,A, A 
particle production for the same collision at the same energy (Table- 2). The cases of k + and k~ 
production in gold-gold reaction at the same energy and at different centralities are reproduced by 
power laws in Fig. 3a and Fig.3b(Table-3). In Fig. 4a and Fig. 4b the production of A and A are shown 
in the same collision at the same energy (Table-4). In Fig. 5a, Fig. 5b and Fig. 5c the production of 
cascade, cascade-bar and omega particles production are shown at different centralities and at the 
same energy (Table-5). The cases of production of neutral kaons and lamda particles in Copper- 
Copper collision at ^/s=200 GeV are plotted in Fig. 6a and Fig. 6b with reckoning of the parameters 
presented in Table-6 and Table-7. 



In Table-8 the values of average transverse momenta (< px >) for different produced secondaries 
in proton-proton and gold-gold collisions have been computed. All these values tally with the similar 
ranges arrived at by experimental measurements. This helps us to obtain for us a consistency 
check-up of the parameter values used for getting fits to the data on < px >-spectra. In Fig. 7a 
and Fig. 7b we see, the lamda-bar to lamda and cascade-minus to cascade-plus particle production 
cross-section ratio as a function of transverse momentum respectively and the ratio gradually fall 
off with increasing values of px- In Fig. 8a and Fig. 8b the nuclear modification factors (Rcp) are 
plotted against transverse momentum for the production of neutral meson and lamda particles in 
copper-copper collision. With the increasing px, the Rcp-v&lues fall off. In addition, the data 
show the Rcp for baryons exhibits a lower fall-off compared with that of mesons in intermediate 
transverse momentum region. The experimental data show that the baryon-meson difference of 
Rcp disappears at higher px- 

The data on production of strange particles described here pertain, in the main, to the 'hard' 
sector of high energy interactions. And it is well-known that Hagedornian power law forms which 
have their roots in the physics of quantum chromodynamics(QCD) describe hard particle production 
in a modestly successful manner for, at least, the light hadrons of which strange K-mesons constitute 
a part. But some of the strange particles have moderately high masses, for which our objective here 
was primarily to check whether this generalized power law form could address the issues of invariant 
py-spectra and some other related observables in a satisfactory manner for all the strange hadrons. 
And the outcome is: this study is strongly affirmative by all indications and yardsticks of actual 
performances. 

The data on strange particles are in general quite sparse. The errors in measurements are also 
in most cases quite considerable. If these limitations are taken into account, the importance of 
this comprehensive work, though entirely phenomenological, assumes some degree of importance. 
Barring these generalized comments on an overall basis, we have some specific observations as well, 
which are also quite well-merited and are being presented hereafter : (i) The power indices for all 
the varieties of strange hadrons lie in the most cases, in the range of 10-12. This is in accord with 



the prescription on the limit set by Brodsky[17j, with the q- values bordering on the values, 2-3. (ii) 
The x 2 / n df values for production of cascade and omega particles have suffered quantitatively due 
to very small values of the number of degrees of freedom, (iii) Thirdly, the q and n values do not 
exhibit any clear centrality-dependence or the mass-dependence of the observed heavy secondaries, 
(iv) However, we cannot ascertain at this moment the properties of these parameters with regard 
to the nature of their energy-dependence(s), if any. (v) The average momentum values of these 
measured heavy strange baryons are found to be quantitatively compatible with other non-strange 
light hadrons, though the expression for the average transverse momentum is not very rigorously 
derived, for which reliability of expression (3) is certainly limited, (vi) Some ratio-values shown by 
Fig. 7a and Fig. 7b are modestly well-described, (vii) However, the values of the nuclear modification 
factor, denoted by Rcp, are not reproduced satisfactorily, especially on the lower-side of the pr- 
values. But this is no wonder, as the used power laws are suited to high-py values as was remarked 
above very concretely, (viii) Still, with one of the simplest approaches, that we have succeeded 
in explaining the characteristics of a large bulk of data on these rare hadrons is certainly quite 
stimulating to and encouraging for us. 

4 Concluding Remarks 

Let us now sum up; 

(i) The used power laws explain quite well the measured data on the observables like, pT-spectra, 
some ratio-behaviours and the nuclear modification factors; so none of the questions related to 
suppression or enhancement is consequential, (ii) The centrality- dependences of the pr-spectra 
of strange hadrons too are well-reproduced, (iii) The essential physical content of the power-law 
models is the modest observance of the p^-scaling (as is reflected in the ~ pt/po term). And in 
terms of the functional efficacy, this model seems, so far, to be at par, if not better, with all other 
numerous existing approaches within the frameworks of either the Hadronic transport models or 
the statistical models[18j-|22j. (iv) The model obviously bears no relationship with the concept of 
the quark-gluon plasma (QGP), which is, still, just a conjecture with so far no clear and concrete 
experimental support, (v) In a conclusive vein, this has to be asserted that we observe nothing too 



strange about 'strangeness' production in high energy interactions. 

The latter two points in the preceding paragraph highlight, in the main, the novelty of this study 
from a global viewpoint against the background of the current trends and streaming ideas in the 
present day Particle Physics. 
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Table 1: Numerical values of the fit parameters of power law equation for keon and lamda 
production in p-p collisions at ^/sjVAr=200 GeV, px = to 5GeV/c 



Sesondaries 


A 


q 


n 


ndf 


K« 


0.563 ±0.023 


3. 108 ±0.025 


15.005 ±0.032 


20.707/17 


K+ 


1.067 ±0.032 


1.581 ±0.055 


10.000 ±0.209 


10.943/9 


K~ 


0.066 ±0.008 


2.895 ±0.026 


15.116 ±0.405 


0.584/6 


A 


0.273 ±0.036 


3.092 ±0.068 


15.007 ±0.074 


18.393/15 


Abar 


0.029 ±0.001 


3.010 ±0.068 


15.016 ±0.068 


15.673/15 



Table 2: Numerical values of the fit parameters of power law equation for lamda and cascade 
particle production in p-p collisions at i/satat=200 GeV, px = to 5GeV/c 



Sesondaries 


A 


q 


n 


ndf 


A 


1.848 ±0.094 


1.306 ±0.061 


10.417± 0.244 


26.883/12 


Abar 


0.541 ±0.026 


1.910 ±0.068 


12.005 ±0.061 


36.675/13 


H - 


0.022 ±0.001 


1.804 ±0.023 


10.033 ±0.119 


21.4/8 


R + bar 


0.063 ±0.002 


2.301 ±0.015 


9.951 ±0.466 


24.246/9 



Table 3: Numerical values of the fit parameters of Power Law equation for keon production (k + , k 
in Au-Au collisions at -y/sArAr=200GeV at different Centrality, pr = to 2 GeV/c 



Secondaries 


Ccntrality(%) 


A 


q 


n 


xL 

ndf 


k+ 


0-5% 


120.114 ± 8.341 


2.508 ±0.057 


10.376 ±0.172 


4.555/14 




20 - 30% 


47.371 ±2.28 


1.998 ±0.029 


9.989 ±0.112 


2.477/15 




40 - 50% 


16.727 ±2.703 


1.843 ±0.082 


10.006 ±0.111 


1.294/12 




60 - 70% 


3468.89 ±253.7 


2.030 ±0.033 


11.981 ±0.154 


3.022/13 


k~ 


0-5% 


435. 675 ±15.85 


2.286 ±0.076 


12.00 ±0.264 


1.032/15 




20 - 30% 


68.462 ±6.505 


2.379 ±0.061 


12.003 ±0.07 


0.246/12 




40 - 50% 


25.454 ±0.901 


2.239 ±0.084 


12.052 ±0.297 


1.566/14 




60 - 70% 


5.058 ±0.283 


2.016 ±0.101 


11.819 ±0.382 


3.356/14 
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Table 4: Numerical values of the fit parameters in Power Law form for (A and A) production in 
Au-Au collisions at y/sNN=200GeV at pt = 1 to 5 GeV/c ranges and for various Centrality-values 
as given below. 



Secondaries 


Ccntrality(%) 


A 


q 


n 


ndf 


A 


0-5% 


63.570 ±4.841 


3.542 ±0.105 


16.124 ±0.140 


10.064/18 




10 - 20% 


772.323 ±48.46 


2.000 ±0.087 


13.870 ±0.259 


5.778/14 




20 - 40% 


141.42 ±11.33 


2.002 ±0.023 


12.815 ±0.099 


7.512/12 




40 - 60% 


34.334 ±2.317 


2.025 ±0.109 


12.820 ±0.308 


7.238/12 




60 - 80% 


19.811 ±0.968 


2.000 ±0.014 


11.869 ±0.235 


14.913/12 


Abar 


0-5% 


198.233 ±0.269 


2.000 ±0.014 


10.765 ±0.269 


7.588/10 




10 - 20% 


29.161 ±1.727 


2.002 ±0.121 


12.003 ±0.072 


17.298/10 




20 - 40% 


122.194 ±4.406 


2.899 ±0.014 


15.001 ±0.054 


3.589/8 




40 - 60% 


25.183 ±0.646 


2.500 ±0.008 


12.888 ±0.028 


9.902/12 




60 - 80% 


10.994 ±0.788 


2.001 ±0.022 


11.995 ±0.088 


31.900/13 



Table 5: Numerical values of the fit parameters in Power Law equation for cascade-minus(H _ ), 
cascade-plus bar(H + 6ar) and Q~ + il + particle production in Au-Au collisions at ^/sATjv = 200GeV 
at different py-values = 1 to 5 GeV/c for various Centrality values. 



Secondaries 


Centrality(%) 


A 


q 


n 


xL 

ndf 


s~ 


0-5% 


213. 632 ±16.06 


1.344 ±0.015 


10.042 ±0.271 


14.593/8 




10 - 20% 


1019.24 ±81.73 


1.001 ±0.010 


10.017 ±0.058 


20.800/6 




20 - 40% 


2978.87 ±155.4 


1.000 ±0.006 


11.003 ±0.001 


6.138/6 




40 - 60% 


16.505 ±0.288 


1.500 ±0.003 


10.019 ±0.055 


0.675/7 




60 - 80% 


3.065 ±0.339 


1.503 ±0.027 


10.001 ±0.111 


16.968/5 


'E.+bar 


0-5% 


1531.25 ±144 


1.001 ±0.012 


9.996 ±0.070 


18.845/7 




10 - 20% 


5248.25 ±468.3 


0.801 ±0.008 


10.011 ±0.053 


22.498/6 




20 - 40% 


735. 176 ±31.68 


0.999 ±0.005 


9.996 ±0.031 


24.258/8 




40 - 60% 


233.981 ±15.39 


0.999 ±0.008 


9.997 ±0.047 


16.173/7 




60 - 80% 


2.993 ±0.188 


1.496 ±0.014 


10.006 ±0.057 


4.049/4 


q- + n+ 


0-5% 


1.248 ±0.094 


3.018 ±0.055 


9.988 ±0.139 


2.073/3 




20 - 40% 


2.970 ±0.229 


1.997 ±0.028 


9.978 ±0.100 


1.304/2 




40 - 60% 


0.342 ±0.037 


2.001 ±0.043 


9.235 ±0.131 


4.755/3 



Table 6: Numerical values of the fit parameters of Power Law equation for keon production(fcg) 
in Cu-Cu collisions at - v /sTVAr=200GeV at different pT-values = 1 to 9 GeV/c 

and for several Centrality domains. 



Secondaries 


Centrality(%) 


A 


q 


n 


ndf 


K 


0-10% 


112.079 ±4.169 


2.020 ±0.035 


12.229 ±0.094 


6.675/21 




10 - 20% 


67.284 ±2.154 


1.921 ±0.031 


11.712 ±0.477 


6.068/21 




20 - 30% 


24.860 ±0.568 


2.879 ±0.030 


14.155 ±0.262 


3.667/21 




30 - 40% 


17.099 ±0.286 


2.565 ±0.070 


13.181 ±0.0.245 


6.230/21 




40-50% 


18.504 ±0.598 


1.985 ±0.036 


11.583 ±0.476 


7.738/21 




50 - 60% 


12.159 ±0.536 


2.000 ±0.006 


11.551 ±0.054 


7.842/21 
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Table 7: Numerical values of the fit parameters of Power Law equation for lamda particle (A) 
production in Cu-Cu collisions at ^/sArjv=200GeV at different Centrality, px = 1 to 7 GeV/c 



Secondaries 


Centrality(%) 


A 


q 


n 


ndf 


A 


0-10% 


1091.23 ±47.1 


2.000 ±0.009 


14.428 ±0.166 


2.635/12 




10 - 20% 


139.73 ±12.91 


2.034 ±0.103 


13.033 ±0.275 


15.067/14 




20 - 30% 


105.009 ±5.974 


2.024 ±0.073 


12.996 ±0.193 


16.425/17 




30 - 40% 


46.450 ±1.523 


2.026 ±0.081 


12.629 ±0.216 


24.731/17 




40-50% 


27.007 ±1.798 


2.001 ±0.018 


12.340 ±0.215 


16.803/15 




50 - 60% 


24.819 ±2.077 


2.029 ±0.023 


13.042 ±0.097 


17.437/14 



Table 8: Calculated values of average transverse momentum < px > for y cm < 0.5 at 
V^j\w=200GeV . 



Nature of Collisions 


Secondaries 


Value of < pt > 


in GeV/c 


p-p Collisions 




0.52 






k+ 


0.45 






k~ 


0.48 






A 


0.52 






A 


0.50 






" — 


0.51 






S + 6ar 


0.86 




Au-Au Collisions 


k+ 


0.68 




(Central) 


k~ 


0.50 






A 


0.54 






A 


0.52 






, ^— 


0.38 






H + 6ar 


0.28 






n- + n + 


0.86 
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Figure 1: Transverse momentum spectra for production of keon (K°, K + , K~), lamda(A), 
lamdabar(A), cascade minus(H~), cascade plus bav(E + bar) particles in pp collisions at \/s = 
200GeV. The experimental data are taken from Ref.[5j. The solid curves are fits for power-law 
model. 
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Figure 2: Transverse momentum spectra for production of cascade(H _ , H + 6ar), lamda and lamd- 
abar particles in pp collisions at - v /s=200GeV. The experimental data are taken from Ref.[23j. The 
solid curves are fits for powerdaw model. 
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Figure 3: Transverse momentum spectra for production of kaon(K + , K~)&t different centrality at 
^s=200GeV in Au-Au collisions. The experimental data are taken from Ref. [24| The solid curves 
are fits for power-law model. 
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Figure 4: Transverse momentum spectra for production of lamda(A) and lamda bar(A) particles at 
different centrality in Au-Au collisions. The experimental data are taken from Ref . }24j The solid 
curves are fits for power-law model. 
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Figure 5: Transverse momentum spectra for production of cascade minus(H~), cascade plus 
bar(H + 6ar) and omega particles(Sl~ + $7 + ) at different centrality in Au-Au collisions at -y/s=200GeV 
. The experimental data are taken from Ref . |25j The solid curves are fits for power-law model. 
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Figure 6: Transverse momentum spectra for production of neutral keons(A^) and lamda(A) particles 
at different centrality in Au-Au collisions at ^/s=200GeV. The experimental data are taken from 
Ref.|24j The solid curves are fits for power-law model. 
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Figure 7: Transverse momentum-dependence spectra of Abar/A and H~ fE, + bar for pp collision at 
- v /s=200GeVThe data type-points are taken from Rcf.|5j. The solid curves or lines are drawn on 
the basis of Power Law Model. 
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Figure 8: Plots of the nuclear modification factor (R cp ) versus pt spectra for Cu-Cu collisions at 
y/s=200GeV. The experimental data are taken from Ref.[24j The solid curves indicate the power- 
law-based description of the data. 
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